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Abstract. Wc prove the existence of a smooth family of non-compact domains ilg C 
j^n+i , n > 1 . bifurcating from the straight cylinder B" x ]R for which the first eigenfunction 
of the Laplacian with Dirichlet boundary condition also has constant Neumann data at 
the boundary: For each s G (— e,e), the overdetermined system 

{Au + Xu = m rig 
u = on dflg 

(Vw, u) = const on dflg 

has a bounded positive solution. The domains fig are rotationally symmetric and periodic 
with respect to the R-axis of the cylinder; they are of the form 

= e R" X R I 11x11 < 1 + scos (y*) +o{s^)^ 

where Ts = Tq + 0(s) and Tq is a positive real number depending on n. For n > 2 
these domains provide a smooth family of counter-examples to a conjecture of Berestycki, 
Caffarelli and Nirenberg. We also give rather precise upper and lower bounds for the 
bifurcation period Tq. This work improves a recent result of the second author. 



1. Introduction and main results 

1.1. The problem. Let f2 be a bounded domain in with smooth boundary, and con- 
sider the Dirichlet problem 




Au + Xu = inf2 
u = on dfl. 



Denote by Ai(f2) the smallest positive constant A for which this system has a solution 
(i.e. Ai(r2) is the first eigenvalue of the Laplacian on Q with Dirichlet boundary condition). 
By the Krein-Rutman theorem, the corresponding solution u (i.e. the first eigenfunction 
of the Laplacian on Q with Dirichlet boundary condition) is positive on Q, and u is the 
only eigenfunction with constant sign in Q, see [TT| Theorem 1.2.5]. By the Faber-Krahn 
inequality, 

(2) Ai(fi) > \i{B'\n)) 
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where B^{Q) is the round ball in R" with the same volume as Q. Moreover, equality holds 
in ([2]) if and only if f2 = B"'{Q), see [S] and [2]. In other words, round balls are minimizers 
for Ai among domains of the same volume. This result can also be obtained by reasoning 
as follows. Consider the functional Q — ?■ Ai(f2) for all smooth bounded domains Q in 
of the same volume, say Vol(f2) = a. A classical result due to Garabedian and Schiffer 
asserts that f2 is a critical point for Ai (among domains of volume a) if and only if the first 
eigenfunction of the Laplacian in fl with Dirichlet boundary condition has also constant 
Neumann data at the boundary, see [9]. In this case, we say that Q is an extremal domain 
for the first eigenvalue of the Laplacian, or simply an extremal domain. Extremal domains 
are then characterized as the domains for which the over- determined system 

{An + An = infi 
u = o on on 

(Vm, u) = const on d^l 

has a positive solution (here u is the outward unit normal vector field along dfl). By a 
classical result due to J. Serrin the only domains for which the system ([3]) has a positive 
solution are round balls, see [18]. One then checks that round balls are minimizers. 

For domains with infinite volume, at first sight one cannot ask for "a domain that 
minimizes Ai". Indeed, with cQ = {cz | 2 G fi} we have 

Ai(cfi) = c~^\i{n), c> 0. 

On the other hand, system ([3]) can be studied also for unbounded domains. Therefore, 
it is natural to determine all domains Q for which has a positive solution. This is 
an open problem. We will continue to call such a domain an extremal domain. In the 
non-compact case, this definition does not have a geometric meaning, except for domains 
which along each coordinate direction of are bounded or periodic. In the case of periodic 
directions, one obtains extremal domains for the first eigenvalue of the Laplacian in flat 
tori, cf. Remark 1 1 . 3 1 b elow . 

Berestycki, Caffarelli and Nirenberg conjectured in [1] that if / is a Lipschitz function on 
a domain f2 in R" such that R'^\fi is connected, then the existence of a bounded positive 
solution to the more general system 

{Am + f{u) = in n 
u = o on on 

(Vu, i^) = const on d^l 

implies that f2 is a ball, or a half-space, or the complement of a ball, or a generalized 
cylinder 5^ x R"~'^' where B'' is a round ball in R'^. In [2U], the second author constructed 
a counter-example to this conjecture by showing that the cylinder x R C R""*"^ (for which 
it is easy to find a bounded positive solution to ([3])) can be perturbed to an unbounded 
domain whose boundary is a periodic hypersurface of revolution with respect to the R- 
axis and such that ([3]) has a bounded positive solution. More precisely, for each n > 2 
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there exists a positive number T^, = T^{n), a sequence of positive numbers Tj — )■ T^, and a 
sequence of non- constant Tj -periodic functions Vj G C^'°(]R) of mean zero (over the period) 
that converges to in C^'"(]R) such that the domains 

= {{x,t) e X R I ||a;|| < 1 + Vj{t)} 

have a positive solution Uj G C'^'°'{Qj) to the problem (Q. The solution Uj is Tj-periodic 
in t and hence bounded. 

1.2. Main results. The goal of this paper is to show that these domains VLj (introduced in 
[20] by the second author) belong to a smooth bifurcating family of domains, to determine 
their approximate shape for small bifurcation values, and to determine the bifurcation 
values T^,(n). Our main result is the following. 

Theorem 1.1. Lei o(R/27rZ) be the space of even 27i-periodic C^'" functions of mean 
zero. For each n > 1 there exists a positive number = T^{n) and a smooth map 

(-£,£) ^ C,^:,,o(W2vrZ) X R 
s ^ iws,Ts) 

with Wq = 0, Tq = T^, and such that for each s G {—e,e) the system ([3]) has a positive 
solution Us G C'^'°'{Qs) on the modified cylinder 

(5) fi, = G R" X R I ||x|| < l + scos(^t ) +5^, ( 



The solution Ug is Ts-periodic in t and hence bounded. 



For n = 2 and for |s| small enough, the bifurcating domains Vis look as in Figured] For 
a figure for n = 1 see Section [8] 



a; G R" 
i 




Figure 1. A domain VLg. 



Notice that for n = 1, the domains Vts do not provide counter-examples to the conjecture 
of Berestycki, Caffarelli and Nirenberg, because R^ \ VLg is not connected. 
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Remark 1.2. From the extremal domains Qs C H^~^^ and the solutions Ug from Theo- 
rem 11.11 we obtain other extremal domains by adding an R'^-factor: For each k > 1 the 
domains fl^ := fig x K-^ are extremal domains in 11"+-'^"'"'^ with solutions t, y) := Us{x^ t) 
(where y G R'^). For instance, in we then have the "wavy cylinder" in Figure [H and 
the "wavy board" obtained by taking the product of the wavy band in Figure |2] with R. 
Notice that R""'"^^'^ \ f2j is connected if and only if n > 2. 

Remark 1.3. The characterization of extremal domains described in Section fLTl more gen- 
erally holds for domains in Riemannian manifolds: Given a Riemannian manifold (M, g), a 
domain C M of given finite volume is a critical point of — )■ Ai(f2), where Ai(f2) is the 
first eigenvalue of the Laplace-Beltrami operator — A^, if and only if the over-determined 
system 

AgM + AM = in Vt 
(6) { u = Q on dn 

g{Vu,u) = const on d^l 

has a positive solution (here u is the outward unit normal vector to dQ with respect to g), 
see [7] and [16]. Theorem 11.11 thus implies that the full tori 



= < {x, t) G R" X R/T, Z 



, „ , /27r \ /27r 

|a;|| < 1 + scos I — t 1 + svsi—t 

are extremal domains in the manifold R" x R/T^Z with the metric induced by the Eu- 
clidean metric. O 



Open problem 1. Are the extremal domains fig K-"" x R/T'sZ (local) minima for the 
functional Vt — )• Ai(f2) ? 

It follows from our proof of Theorem 11.11 and from the Implicit Function Theorem that 
the family is unique among those smooth families of extremal domains bifurcating from 
the straight cylinder that are rotationally symmetric with respect to R" and periodic with 
respect to R. A much stronger uniqueness property should hold. Indeed, the existence 
problem of extremal domains near the solid cylinder, say in R^, is tightly related to the 
existence problem of positive constant mean curvature surfaces near the cylinder, see Sec- 
tions [2] and [3l Any positive constant mean curvature surface with two ends (that is properly 
embedded and complete) must be a Delaunay surface, by a result of Korevaar, Kusner, 
and Solomon, [13]. We thus ask: 

Open problem 2. Assume that Vt is an unbounded extremal domain in R"^-^ that is 
contained in a solid cylinder. Is it then true that Vt belongs to the family VLg ? 

We also determine the bifurcation values = T^{n). It has been proved in [20] that 
T*{n) < In particular, T*(r2) — > as n — oo. We shall show in Section [8] that 

T,(l) = 4. Fix now n > 2 and define v=^. Write for T,(n). 



BIFURCATING EXTREMAL DOMAINS FOR THE FIRST EIGENVALUE OF THE LAPLACIAN 5 



Theorem 1.4. Let Jy\ (0, +00) — )■ ]R &e the Bessel function of the first kind. Let jy he its 
smallest positive zero. Then the function sJy_i{s) + Jy{s) has a unique zero on the interval 
(0, jj,), say p^, and 

27r 

T^, = 

- pI 

In particular, 

Furthermore, the sequence Ty is strictly decreasing to 0. 

The numbers Ty for z/ < 10 are given in Section O In particular, for n = 2, 3 and 4 
(corresponding to the bifurcation of the straight cylinder in R^, and R^) the values of 
Ty are 

To ^ 3.06362, Ti ^ 2.61931, Ti ^ 2.34104. 

Open problem 3. Is the bifurcation at T^{n) sub-critical, critical, or super- critical? In 
other words, 9,(T,)|,^o < 0, 9,(T,)|,^o = 0, or 9,(T,)|,^o > ? 

The paper is organized as follows. In Section |2] we show how the existence of Delaunay 
surfaces (i.e., constant mean curvature surfaces of revolution in that are different from 
the cylinder) can be proved by means of a bifurcation theorem due to Crandall and Rabi- 
nowitz. We will follow the same line of arguments to prove Theorem II. II in Sections |3] to [HI 
In Section [9] we prove Theorem 11.41 on the bifurcation values T*(n). 

Acknowledgments. Most of this paper was written in June 2010, when the second 
author visited Universite de Neuchatel. The second author is grateful to Bruno Col- 
bois and Alexandre Girouard for their warm hospitality. The first author thanks FRU- 
MAM for its hospitality during the workshop "Problemes aux valeurs propres et problemes 
surdetermines" at Marseille in December 2010. We both thank Frank Pacard for helpful 
discussions and for kindly allowing us to include the exposition in Section [2l 

2. The Delaunay surface via the Crandall-Rabinowitz Theorem 

Our proof of Theorem 1 1.1 1 is motivated by the following argument that proves the existence 
of Delaunay surfaces by means of the Crandall-Rabinowitz bifurcation theorem. The 
material of this section was explained by Frank Pacard to the second author when he was 
his PhD student. 

We start with some generalities. Let S be an embedded hypersurface in R""*"^ of codi- 
mension 1. We denote by // its second fundamental form defined by 

II{X,Y) = -{VxN,Y) 

for all vector fields X, Y in the tangent bundle T S. Here N is the unit normal vector field 
on S, and (■, ■) denotes the standard scalar product of R"""*"^. The mean curvature i7 of S 
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is defined to be the average of the principal curvatures, i.e. of the eigenvalues /ci, . . . , /c„ of 
the shape operator A : T S — y T S given by the endomorphism 

{AX,Y) = -IIiX,Y). 

Hence 

n 
i=l 

Given a sufficiently smooth function w defined on S we can define the normal graph 
of w over E, 

= {p + w{p) N{p) G I P G S} , 

and consider the operator w i— )■ HiY^uy) that associates to w the mean curvature of E^. 
The linearization of this operator at w = is given by the Jacobi operator: 

where g the metric induced on E by the Euclidean metric and — is the Laplace-Beltrami 
operator on E. All these facts are well-known, and we refer to [2] for further details. 

In 1841, C. Delaunay discovered a beautiful one-parameter family of complete, embed- 
ded, non-compact surfaces D^^ in R^, a > 0, whose mean curvature is constant, see 
These surfaces are invariant under rotation about an axis and periodic in the direction of 
this axis. The Delaunay surface D^r can be parametrized by 

X,{9,t) = {y{t) cos 9, y{t) sine, z{t)) 

for {9, t) & X R, where the function y is the smooth solution of 



and z is the solution (up to a constant) of 



At) '■^'<" + " 



When cr = 1, the Delaunay surface is nothing but the cylinder Di = S*^ x R. It is easy to 
compute the mean curvature of the family D^j and to check that it is equal to 1 for all a. 
One can obtain each Delaunay surface by taking the surface of revolution generated by 
the roulette of an ellipse, i.e. the trace of a focus of an ellipse i as ^ rolls along a straight 
line in the plane. In particular, these surfaces are periodic in the direction of the axis of 
revolution. When the ellipse I degenerates to a circle, the roulette of I becomes a straight 
line and generates the straight cylinder, and when a — >■ 0, D^j tends to the singular surface 
which is the union of infinitely many spheres of radius 1/2 centred at the points (0,0, n), 
n G For further details about this geometric description of Delaunay surfaces we refer 

to [g. 
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We now prove the existence of Delaunay surfaces by a bifurcation argument, using a 
bifurcation theorem due to M. Crandall and P. Rabinowitz. Their theorem apphes to 
Delaunay surfaces in a simple way. We shall use the same method to prove Theorem II. 1[ 
The phenomenon underlying our existence proof of Delaunay surfaces is the Plateau- 
Rayleigh instability of the cylinder, [T7] . 

Consider the straight cylinder of radius 1, in cylindrical coordinates: 

Ci = {(p, 6, t) e (0, +oo) X 5^ X R I p = 1}. 
Let w be a C^-function on x R/27rZ. In Fourier series, 

w{9, t) = (aj cos(j 6) + Pj sin(j 6')) (a^ cos [k t) + bk sin {k t)^ . 

j,k>0 

If w{6,t) > —1 for all 6,t, we consider, for each T > 0, the normal graph Cj^^ over the 
cylinder Ci of w rescaled to period T, 

Cl+^ := !^{p,e,t)e{0,+oo)xS'xR I p = l + nj(9,'^?jy 

Define the operator 

where H is the mean curvature. Then F{w, T) is a function on S*^ x R of period T in the 
second variable. Therefore, 

(7) Fiw,T){9,t) := F{w,T) [o.^t 

is a function on x R/27rZ. Note that F(0,T) = for all T > 0, because for w = 
the surface is the cylinder Ci whose mean curvature is 1. If we found a non-trivial 

solution {w,T) of the equation F{w,T) = 0, we would obtain a constant mean curvature 
surface different from Ci. In order to solve this equation, we consider the linearization of 
the operator F with respect to w and computed at {w,T) = (0,T). As mentioned above, 
the linearization of the mean curvature operator for normal graphs over a given surface 
with respect to w computed at = is the Jacobi operator. Since the Laplace-Beltrami 
operator on Ci (with the metric induced by the Euclidean metric) is —dg — d^, and since 
the principal curvatures ki of Ci are equal to and 1, we find that 

D^F(0,T) = • 

For each j, /c G N U {0} and each T > 0, the four 1-dimensional spaces generated by the 
functions 




cos(j^) cos{kt), cos{j6) sm{kt), sin{j6) cos{kt), sin{j6) sin{kt) 
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are eigenspaces of F(0, T) with eigenvalue 




Clearly, 

• cTj, fe(T) ^ for all T > if j > 2, or if j = 1 and /c > 1; 

• (Ti'o(T) = for all T > 0; 

• (7o.fc(T) = only for T = 2TTk and k > 1; moreover o"o.fc(^) changes sign at these 
points. 

It follows that Ker DwF{0,T) is 2-dimensional (spanned by cos^, sin^) if T > and 
T ^ 27rN, and that Ker Dy^F{0,T) is 4-dimensional (spanned by cos^, sin^, cos(/ct), 
sin(H)) if T e 27rN. 

We will now bring into play an abstract bifurcation theorem, which is due to Crandall 
and Rabinowitz. For the proof and for many other applications we refer to |T2l [19] and to 
the original exposition j3]. 

Theorem 2.1. (Crandall— Rabinowitz Bifurcation Theorem) Let X and Y be Ba- 

nach spaces, and let U G X and A C R 6e open subsets, where we assume G t/. Denote 
the elements of U by w and the elements of A by T. Let F: U x A ^ Y be a C°° -smooth 
function such that 

i) F(0,T) = for allT eA; 

ii) Ker F{0, Tq) = Hwo for some Tq G A and some wq E X \ {0}; 

iii) codimlm F(0, To) = 1; 

iv) DtD^F{0,To){wo) i lmD^F{0,To). 

Choose a linear subspace X G X such that HwqQX = X . Then there exists a -smooth 
curve 

(-£, e) ^ X X R, {w{s),T{s)) 

such that 

1) w{0) = and T{0) = Tq; 

2) s {wo + w{s)) e U and T{s) e R; 

3) F{s{wo + w{s)),T{s))=0. 

Moreover, there is a neighbourhood Af of (0, Tq) E X xH such that {s (wq + w{s)) , T(s)) } 
is the only branch in M that bifurcates from {(0,T) | T G A}. 

The theorem is useful for finding non-trivial solution of an equation F[x, A) = 0, where x 
belongs to a Banach space and A is a real number. It says that under the given hypothesis, 
there is a smooth bifurcation into the direction of the kernel of Dy^F for the solution of 
F(x, A) = 0, and that there is no other nearby bifurcation. 

In order to apply Theorem 12.11 we now restrict the operator F defined in ([7j) to functions 
that are independent of 9 (so as to get rid of the functions cos 6', sin 6' in the kernel of 



BIFURCATING EXTREMAL DOMAINS FOR THE FIRST EIGENVALUE OF THE LAPLACIAN 9 



D^F{0, T)) and that are even (so as to have a 1-dimensional kernel for T G 27rN). We can 
also assume that the functions w have zero mean. In other words, we look for new constant 
mean curvature surfaces among deformations of Ci that are surfaces of revolution, even in 
the t-direction. We hence consider the Banach space 

of even 27r-periodic functions of zero mean whose second derivative is Holder continuous. 
Moreover, define the open subset U = {w & X \ w(t) > —1 for all t} of X, and the Banach 
space 

Y = CXo(W2vrZ). 

Furthermore, chose A = (0, +oo) C R. Then the operator F defined as above restricts to 
the operator 

F:UxA-^Y. 

With 

crfc(T) := (To,fc(T) = 




its linearization with respect to w at Tq := 27r is 



Dy,F{0,To) f^flfc cos(A;t) j = ^(Tfc(To)afc cos(/ct) = ^(A;^ 

\fc>l / k>l k>l 



1) ttk cos(fct). 



Hence, 

KerL)^F(0,To) = R cost. 

Moreover, the image ImDm F{0, Tq) is the closure of ®fc>2 Rcos{kt) in Y; its complement 
in Y is the 1-dimensional space spanned by cos t. Finally, 



dT 



cost = cost ^ lmDwF{0,To) 

T=To ^TT 



With wq = cost and X the closure of ^^>2 ^'^'^^(/ct) in X, the Crandall-Rabinowitz 
bifurcation theorem applies and yields the existence of C°°-smooth curve 

{-e,e) ^ X xR, s ^ {w{s),T{s)) 

such that 

1) w{0) = and r(0) = Tq; 

2) F{s{wo + w{s)),T{s))=0, 

i.e. (by the definition of the operator F) the existence of a C°°-smooth family of surfaces 
of revolution that have mean curvature constant and equal to 1, bifurcating from the 
cylinder Ci. That these surfaces are Delaunay surfaces follows from Sturm's variational 
characterization of constant mean curvature surfaces of revolution, [H [5] . 
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Remark 2.2. The boundaries of the new domains Qg C described in Theorem 11.11 
are not Delaunay surfaces (at least not for |s| small). Indeed, Delaunay surfaces bifurcate 
from the cylinder at Tq = 27r, while the domains Qg bifurcate from the cylinder at T*(2) ^ 
3.06362. 

3. Rephrasing the problem for extremal domains 

We want to follow the proof of the existence of Delaunay surfaces given in the previous 
section in order to prove the existence of a smooth family of normal graphs over the 
straight cylinder such that the first eigenfunction of the Dirichlet Laplacian has constant 
Neumann data. In this section we recall the set-up from [20], where the second author 
studied the Dirichlet-to-Neumann operator that associates to a periodic function v the 
normal derivative of the first eigenfunction of the domain defined by the normal graph of v 
over the straight cylinder, and computed the linearization of this operator. The novelty of 
this paper is the analysis of the kernel of the linearized operator; it will be carried out in 
Sections m to [71 

The manifold ]R/27rZ will always be considered with the metric induced by the Euclidean 
metric. Motivated by the previous section, we consider the Banach space C^^^j^ Q(M/27rZ) of 

even functions on M/27rZ of mean 0. For each function v G C'^^'^^^o{M./27t'Ij) with v(t) > —1 
for all t, the domain 



Cf^^ := <j {x,t) e R" X R/TZ \ < ||x|| <l + v{—t 



2tt 

Y 

is well-defined for all T > 0. The domain Cj^^ is relatively compact. According to standard 
results on the Dirichlet eigenvalue problem (see [ID]), there exist, for each T > 0, a unique 
positive function 

and a constant A = r G IR such that is a solution to the problem 



(8) 

which is normalized by 



A0 + A0 = in 



(9) / ) dvol 



T 



on dCf^, 



2 



Furthermore, (p and A depend smoothly on v. We denote 0i := </)o,t and Ai := Ao,t- Notice 
that 01 does not depend on the t variable and is radial in the x variable. (Indeed, 0i is 
nothing but the first eigenfunction of the Dirichlet Laplacian over the unit ball B" in 
normalized to have L^-norm ^.) We can thus consider 0i as a function of r := and 
we write 

(10) y;i(r)=0i(x). 
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We define the Dirichlet-to-Neumann operator 

where u denotes the unit normal vector field on dCjj^^ and where (p = 4'v,t is the solution 
of dH]). The function 

F{v,T): (9Cf+,, = d{W) X R/TZ ^ R 
depends only on the variable t G R/TZ, since v has this property. It is an even function; 
indeed, v is even, and hence 0„,t is even, since the first eigenvalue \v,t is simple. Moreover, 
F{v,T) has mean 0. We rescale F and define 

F{v,T){t) = F{v,T) i^t 



2n 

Schauder's estimates imply that F takes values in Cj!^"^ Q(IR/27rZ). With 
U ■= [ve C,^;°„o(W2vrZ) I v{t) > -1 for all t} 

we thus have 

F: t/x (0,+oo) ^ Clf,^^,{R/27rZ). 

Also notice that -F(0, T) = for all T > 0, and that F is smooth. 
The following result is proved in [20] . 

Proposition 3.1. The linearized operator 

Ht := D..F{0,T):Ct:^^,{R/2nZ) ^ Cl:,^,iR/27rZ) 

is a formally self adjoint, first order elliptic operator. It preserves the eigenspaces 

Vk = R cos(H) 

for all k and all T > 0, and we have 

(11) HT{w){t) = (^dr^ + 9,^01 ■ W 

where is the unique solution of 

A^ + Ai^/^ = zn Cf 

i/j = -9,01 ■ w{2TTt/T) on del 

which is L'^{Cf) -orthogonal to (pi, and where r = \\x\\. 
Write 

w{t) = Ofc cos(/ct). 

k>l 

Since Ht preserves the eigenspaces. 



(12) 



HT{w){t) = ^ CTfc (T)afc COS (fct). 



fc>i 
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We use f lTTjl and (fT2|) to describe (JkiT) as the solution of an ordinary differential equation: 
The solution ip of f[T^ is differentiable, and even with respect to x for fixed t. Therefore, 
for each t, the derivative of ip with respect to r vanishes at 0: drip\r=o = 0. Hence, 

(13) CT,{T) = + 

where for n > 2, Ck is the continuous solution on [0, 1] of the ordinary differential equation 

U + ^dr + X.-^-j^j j c, = 

such that Ck{l) = —ip[{l), while for n = 1, Ck is the solution on [0, 1] of the ordinary 
differential equation 




such that Cfc(l) = —(p[{l) and c'f^{0) = 0. Notice that for all /c > 1 and all n > 1 

a,{T) = a, . 

Our next aim is to find an explicit expression for the function ai in order to describe the 
spectrum of the linearized operator, to read off its kernel, and to find the codimension of 
its image. We first consider the case n > 2, for which we need Bessel functions. The case 
= 1 is discussed in Section [HI 



4. Recollection on Bessel functions 

In what follows we shall use several basic properties of Bessel functions. For the readers 
convenience, we recall the definition of the Bessel functions Jr and Ir, and state their 
principal properties. For proofs we refer to [211 Ch. III]. 

4.1. The functions Jr. For r > the Bessel function of the first kind Jr'- R — )■ R is the 
solution of the differential equation 

(14) y'\s) + s y\s) + (s^ - r^) y{s) = 

whose power series expansion is 

m=0 

We read off that 

(16) Jo(0) = 1, Jr{0) = for all r > 0. 

The power series f|T5|) defines a solution Jr'. (0, oo) — R of f|T^ also for r < 0. If r = n is 
an integer, then 

jjs) = (-irjjs) 
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and Jn is bounded near 0. If r is not an integer, then the function Jr{s) is bounded near if 
r > but diverges as s — if r < 0. The functions Jt{s) and J-t{s) are therefore hnearly 
independent, and hence are the two solutions of the differential equation on (0, oo). 
For all r G R and all s > we have the recurrence relations 

(17) Jr-l{s) + Jr+l{s) 

(18) J._l(s) - Jr+l{s) 

(19) sj;(s) + rJ,(s) 

(20) sj;(s) - rJ^{s) 

Another important property that we will use often is that the first eigenvalue Ai of the 
Dirichlet Laplacian on the unit ball of R", n > 2, is equal to the square of the first positive 
zero of Jy for u = Notice that Ai depends on n. Moreover, the function J^, is positive 
on the interval (0, y/Xi), and J^(a/Ai) < 0. 



= y^r(s), 

= 2J:(s), 

= -Sjr+l(s). 



4.2. The functions For r G R the modified Bessel function of the first kind : R — )■ R 
is the solution of the differential equation 

s'y"{s) + sy'{s)-{s^ + T^)y{s) = 

whose power series expansion is 



(21) lAs) = J2 



ml r(r + m + 1) 

m=0 

We read off that It{s) > for all r G R and s > 0, and that 

(22) /o(0) = 1, J^(0) = for all r > 0. 

Comparing coefficients readily shows that for all r G R and all s > we have the recurrence 
relations 

(23) Ir^l{s) - Ir+l{s) = —Ir{s), 

S 

(24) Ir^,{s) + Ir+l{s) = 2/:(s), 

(25) slUs) + Tlr{s) = slr.,{s), 

(26) slUs)-Tlr{s) = Sir+lis). 

We shall also make use of the asymptotics 

(27) lim 4^ = 1. 
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5. A FORMULA FOR cri(T) WHEN n>2 

In this section we begin our analysis of the first eigenvalue (Ji{T) of the linearized oper- 
ator Ht- We assume that n > 2 throughout. To simplify the notation, we denote the 
previously defined function ci by c. Recall that for n > 2, 

a,{T) = c'(l) + ^'/(l) 

where c is the continuous solution on [0, 1] of the ordinary differential equation 

(28) [d^r + ^dr + X,-(^^yy = Q 

such that c(l) = —(p[{l). We shall distinguish three cases, according to whether the term 



T 

is negative, zero or positive. Recall that Ai depends on n. In order to simplify notation, 

n-2 



we put 1/ = and write for Ai = \\irL). As mentioned in the previous section, a/\ 
is the first zero of Jy. Denote 



O'Lcft(r) 


if 


T < 






if 


T = 




crRight(7') 


if 


T > 





and /i = y^. We shall find an explicit expression for o"i(T). For T > denote 



(29) 

5.1. A formula for o"Lcft- Assume that T < fi. This allows us to define 

(30) e = 

We rescale the function c by defining 

CIS) = c \ - 




In view of (!28|) . c is the continuous solution on [0,^] of 



with c{C,) = —ip'i{l). This equation is very similar to a modified Bessel equation. In order 
to obtain exactly a modified Bessel equation, we define the function c by 
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Note that — z/ < because n>2. Hence c is the continuous solution on [0, ^] of 

+ - + c = 



S \ S'^ 



with c{C,) = —^^ ip'i{l). The solution of this ordinary differential equation is given by 
a/j/(s), where the constant a (depending on u and T) is chosen such that 

aJ.(0 = -e>'i(l). 
Returning to the function c, we get 

c{r) = -^^"^^^(^^) 
and from ([13]) and ([29]), using the identities ([M]), ([25]) and ([26]), we obtain 

aLeft(T) = Vi(l)y^^(e4-i(0-2///.(0+a+i(0)+¥'i(l) 

(31) = ^iW-v^'id)^^^- 

To better understand ai^chiT) we shall need the values v^'(l) and (p"{l). From ffTOl) and the 
definition of (pi we have that ipi is the continuous solution on [0, 1] of 

n — 1 \ 
dl + dr + KWi = Q 

such that V5i(l) = 0, with normalization 

1 1 

<y5^(r) dr 



27rVol(5"-i)' 
We rescale the function tpi and define 

Hence, (pi is the continuous solution on [0, ji,] of 
(32) {^l + '^d^ + ^^^ 

with ^pi{ju) = and normalization 

,^2 



(^i(s) ds 







27rVol(5"-i)' 

Equation fl32]) is very similar to a Bessel equation. In order to obtain exactly a Bessel 
equation, we define the function (fi by 
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Since —V < because n > 2, we get that 0i is the continuous solution on [0,ju] of 

01 = 



s \ 



with = and normahzation 



.2-n ,-2/,N J . Ji^ 



27rVol(5"-i) 

The solution of this ordinary differential equation is KnJy{s), where the constant is 
chosen such that 



Jo 



27rVol(S"-i)' 
Returning to the function </)i, we get 

It follows that 

Since Jv{iu) = we obtain 

(33) = «:njr+'J:(j.). 
Furthermore, 

and hence 

cp'lil) = n^j-^+'[-2vJ'M+j,j:{j,)). 
To rewrite this further note that, by flTSl) . 

2J:{s) = Jl_M-JUi{s). 
Together with ([20]) and ([HD we find 

2sJ:{s) = sJl^,{s) - sJl^,{s) 

= ({v - 1) J._i(s) - sJ,(s)) - [-{v + 1) J.+i(s) + sJ,(s) 
At s = we obtain, together with f[T7|) and ffTSl) . 

'^JuJ'Xiv) = Jv+i{ju) - Jv-i{jv) = -2Jl{j^). 

Altogether, 

(34) ^;'(1) = -«:„jr^M2i^ + l)J:(j.). 
In view of (!3T!) . (133|) and (134|) the function (TLcft(r) is equal to 

e/.+i(o 



(35) aLeft(T) = j;-+V:(j.) ((2z/ + 1) + ^ 



(0 



BIFURCATING EXTREMAL DOMAINS FOR THE FIRST EIGENVALUE OF THE LAPLACIAN 17 

Using also ( 125|) and fl26|) we can rewrite this as 

(36) aLeft(T) = -n^j--^^\Jl(j,)^l + ^-My 

Since k^, jy are positive, Jl{ju) is negative, and the functions I^, are positive at all ^ > 0, 
formula ( |36l) implies 

Lemma 5.1. In the interval of definition (0,/i) of the function o"Lcft; have 

CTLeft(T) > 0. 

Moreover, by (|2T|) we have 

Since ^ — ?■ as T /i by f l30p . we find together with that for all > 0, 

ai(/i) = limaLeft(T) = j-+V:(j,) (2z/ + 1 + 2 ^^^^) . 

In particular. 

Lemma 5.2. cri(/i) > 0. 

5.2. A formula for crRight- We follow the reasoning that we used to find a formula for the 
function ai,eit{T). We skip the technical details. Assume that T > fi. This allows us to 
define 



(37) p=^/A.-(^— 
The function c(s) := s'^c(-) is the continuous solution on [0,p] of 



,2 



9^ + -9.-^1 + ^ 







s \ 

with c(p) = —p'^ip'i{l). The solution of this ordinary differential equation is given by 
/3 Ju{s), where the constant (3 (depending on z/ and T) is chosen such that 

/3^.(p) = -P>;(i)- 

Returning to the function c, we get 

c(r) = -^r~^Jy(pr) 
MP) 

and from ( !T3|) and ( 129|) . using the identities ( !T8|) . ( IT9|) and (!20|) . we obtain 

f^Right(T) = -^'i(l)j^^(p-^z^-i(p) - 2z^'//.(p) - pJu+i{p)) + ¥^1(1) 

(38) = ¥^;'(i) + ^;(i)^^^±i^. 
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In view of fl33|) and f p^ this becomes 

(39) an^MT) = -^n3~''^'J'M{{2u+l)-P-^j^ 

where we used the identities f|T9|) and f l20|) to get the second equahty. 

6. Study of the derivative of (yi{T) 
Throughout this section we assume again that n>2. We start with 
Lemma 6.1. The function ai : (0, oo) — )■ R has the asymptotics 

hm cTi(T) = +00 and hm ai{T) = —oo. 

T-s-O r-s>oo 

Proof. The first asymptotics is aheady proven in [20]. We give an easier proof: By (150]) 
we have ^ — )■ oo as T — )■ 0. Using fl5^ and (1271) we therefore find 

hm (7i(T) = hm ^ '^~^}jP = hm £ = oo. 

To prove the second asymptotics, we read off from (1371) that p /• \/\^ = jj, as T — )■ oo. As 
is well-known, j^, < j^+i (see e.g. [211 §15-22]). Therefore J^+i{j^) > 0. Together with ( 139|) 
we thus find 

hm (Ji{i ) = — lim — — - — = — oo. 
as claimed. □ 

It is shown in [201 P- 336] that the function ai is analytic and hence differentiable. For our 
purposes, it would be enough to know that ai has exactly one zero and that cr[{Ty) ^ 0. 
This follows from Lemma 15.11 Lemma 15.21 Lemma 16.11 and Lemma 16.51 below, that states 
that (y'i{T) < for all T G (/i, oo). We shall prove a somewhat stronger statement, namely 
that a[{T) < for all T G (0, oo). 

Proposition 6.2. Let n > 2. The function ai : (0, oo) — j- R has negative derivative. 
Moreover, ai has exactly one zero, say T^. 

Proof. We show that aLcft has negative derivative (Lemma 16.31) . that alight has negative 
derivative (Lemma 16. 5p . and that (j[{fi) < (Lemma 16. 7p . The fact that ai has exactly 
one zero then follows together with Lemma 16.11 

Lemma 6.3. cr'^citC^) <^ '^^^ ^ ^ (0,/i). 



Proof. Recall from (133 p that 



Vi(i)>o. 
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Set f{s) = ^tJ^. In view of we need to show that ^/(^(T)) < for all T G (0,/i). 
Since ^f{^{T)) = /'(^(T)) {'(T) and ^(T) < for all T G (0,/x), this is equivalent to 

(40) f'{s) > for all s G (0,cx)). 

By ( |T9l) we have = — (i^ + + s/;^ and sl^ = —ul^ + sl^^i. Therefore, 

The lemma now follows from the following claim. 

Claim 6.4. J^(s) > I^^i{s)Iu+i{s) for all u e R and all s > 0. 

Proof. In view of (1221) we have /^(O) > (0)1,^+1 (0) for all u > 0. It therefore suffices 
to show that for all s > 0, 

as as 

Multiplying by s, we see that this is true if and only if 

(41) 2I,sIl > sJ;„iJ,+i + h^isll^^. 
In view of ([24]), (EE]), ([25]) we have 

2sll = sl^_i + Sly+l 

sll_i = {y - l)Iu-i + sly 
sll+i = -{y + 1)4+1 + sly. 
Therefore, (j^T!) holds if and only if 

Slu^lly + Sljy + l > {U - + l + Sljy+l - (^ + l)ly^lly+l + Sly^Jy 

i.e., 

> -2/,_iJ,+i 

which is true because Iv{s) > for all G IR and s > 0. □ 

Lemma 6.5. o"^^ight(^) < for all T G (yU, oo). 

Proof. Recall that —(p[{l) > 0. Note that the function 



2vr^^ 



(/i,oo) -> (0, j^), p(T) = yA^-(^Y 

is strictly increasing. Set h{s) = ^^^jj^- In view of ( 138]) we need to show that 
(42) h'{s) > for all s G (0,j^). 

Since jj, is the first positive zero of J^, we see as in the proof of Lemma [6.31 that (l42l) is 
equivalent to 

Claim 6.6. J^{s) > J^_i(s) Ji,+i(s) for all s G {0,j,y). 
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Proof. Let again ju-i, ju, ju+i be the first positive zero of Ju-i, Ju, Ju+i, respectively. 

(2) 

Moreover, denote by the second positive zero of J,,^i. Then 

(43) ju-i < ]u < Ju < ji-i, 

see e.g. [211 §15-22]. It follows from the power series expansion (fTSll that 

(44) J.(s)>0 forse(0,j,). 

Assume first that s G [ju-i,ju)- Then (H3l) and (jSj) show that Ju{s) > 0, Ju^i{s) < 0, 
Ju+i{s) > 0, whence the claim follows. Assume now that s G (0,jjy_i). In view of ([T6|l we 
have J^(0) > Jy_i(0) Jy+i(0) for all u >0. It therefore suffices to show that 

d d 

(45) — J2 > — (J^_iJ^+i) on (0,j^_i). 

as as 

Using ( !T8|) . (!20l) and (ITQ!) we see as in the proof of Claim [63] that (l45l) is equivalent to 

> -2J,„i(s) J,+i(s) 
which is true because J^-i and J^+i are positive on (0, jjy_i). □ 

To complete the proof of Proposition 16.21 we also show 
Lemma 6.7. cr^(/i) < 0. 
Proof. Since the function cti is smooth, 

^'i(/^) = f^Right(^)- 
For T > fi we have cr'^(T) = h'(^p(T)^ p'(T). We compute 

M^) = ;^(T)T^ 

and 

S (J^ — J^_iJj,^i) 



(46) /^'(s) ^2 

Since lim^-x^^ p(^) = we obtain 

a[{fi) = limaRi ht(^) = ^<^i(l) (l-l™-^^^ 



In view of the power series expansion (1151) . 



Therefore, 



Ju^iJu+i r(z/ + 1) 

lim — = ^, < 1 for all u > 



and thus o"'^(/i) < 0. □ 
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7. Extremal domains via the Crandall-Rabinowitz theorem 

We are now in position to prove our main result when n > 2: The hypotheses of the 
Crandall-Rabinowitz bifurcation theorem are satisfied by the operator F defined in Sec- 
tion [31 For n > 2, Theorem 11.11 follows at once from the following proposition and the 
Crandall-Rabinowitz theorem. As before, v = 

Proposition 7.1. Forn > 2, there exists a real number T^{n) = T^, such that the kernel of 
the linearized operator F{0,T^) is 1- dimensional and is spanned by the function cost, 

KerD„F(0,T^) = R cost 

The cokernel of D^F{0, TJ) is also 1- dimensional, and 

DTD,F{0,T,){cost) i ImD,F(0,T,). 

Proof Let v E C'^;:,^,iR/2nZ), 

V = ttk cos{kt). 

k>l 

We know that 

(47) D,F{0,T) = Y,^kiT)ak cos(A;t). 

k>l 

Let Vk be the space spanned by the function cos{kt). By Proposition |621 the function o"i(T) 
has exactly one zero T^. By (147]) . the line Vi belongs to the kernel of F(0, T^). Moreover, 
Vi is the whole kernel, because for k > 2 we have 

(TkiT,) = ^ 

(because Ty is the only zero of ai). By (H7|) and since F{Q,Ty) is elliptic, the image of 
F(0, Ty) is the closure of 

k>2 

in Cg^gj^ Q(]R/27rZ)), and its codimension is equal to 1. More precisely, 

Cg^:,,o(W27rZ)) = lmD,F{0,Ty)(BVi. 

Again by (H71) . 

DtD,F{0,T){v) = Y,^'k{T)ak cos(fct) 

k>l 

and in particular 

DTD^F{0,Ty){cost) = a[{Ty) cost ^ ImL)„F(0,T^) 
because (T[{Ty) < by Proposition 16.21 This completes the proof of the proposition. □ 
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8. The problem in 

Assume that n = 1, i.e., the ambient space of the cyhnder Cf is R^. Recall from Section E] 
that in this case, 



a,iT) = c'(l) + v^'/(l) 



where c is the solution of 
(48) 



c = 0, 



with c(l) = —(p[{l) and c'(0) = 0, where (pi is the first eigenfunction of the Dirichlet 
problem on [—1,1] normalized to have L^-norm (Here and in the sequel, c denotes 
again the function ci.) For ipi and Ai we thus have 



Ai 



TT 



and 



1 /TT 

COS — r 
2tt V2 



Hence 



= Vf and ^;'(1) = 0. 
Lemma 8.1. The only zero of the function cri(T) is at T = A. Moreover cr^(4) < 0. 
Proof We abbreviate a(T) := Ai - {^Y = (f )^ - (^)^ The solution to (gH]) is 



c(r) 



cosh '\y—a{T) r 



8 cosh ^/-"(T) 



if TG (0,4), 
if T = 4, 
if T G (4, oo) 



Hence, 



f V-a(T) tanh v/-«(r) if Tg (0,4), 

if T = 4, 
if T G (4,oo). 



I y a(T) tan y a(T) 



ai(T) = c'(l) 



In particular, cri(T) > on (0,4) and o"i(T) < on (4, oo). It remains to show that 
a[{A) < 0. 

For T > 4 define h{T) := ^a{T). Then 

a[{T) = -./|^(/i(T)tan/i(T)) = h'{T) {tan h{T) + h{T){l + tan\h{T)))) . 
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Since cri(T) is smooth on (0, oo) and since lini^ h{T) = and h'{T) = , we find 

a[{A) = \im^ h'{T) {tan h{T) + h{T)) 

f lim h'(T)2h(T) 
-Jf lim a'{T) = -Jf ^ < 0. 



□ 



Remark 8.2. A computation shows that c'l^T) < for all T G (0, oo). 

Using the previous lemma, the proof of Proposition 17.11 applies also for n = 1, and we 
obtain 



Proposition 8.3. Proposition 7.1 is true also for n = 1 and T*(l) = 4. 

Together with the Crandall-Rabinowitz theorem we now obtain our main Theorem 11.11 
also for n = 1. Figure [2] shows the shape of the new extremal domains in R^. 




Figure 2. A domain fi^ c R^. 



9. Estimates on the bifurcation period 

Recall from Section |8] that T^,(l) = 4. In this section we study the bifurcation values 
Ti, = T^{n) for n>2, and in particular prove Theorem 11.41 

We recall that Jl{ju) 7^ 0, and from fH2|) that the function h{s) = ^'^j'^^j^^^ is strictly in- 
creasing on (0, j;,) from to oo. By fl39|) the unique zero T^, of CTRight is therefore determined 
by 

(49) := 
and 
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In other words, the bifurcation value is 
(50) 



27r 



a/a: 



Pt 



where is the unique zero on (0, jj,) of sJ^+i — {2u + \)Ju or, by f lT7|l . of sJ^^i + J,^. 

For fixed the value p^, and hence can be computed by the computer (using, for 
instance, Mathematica). The first few and some larger values of (rounded to five decimal 
places) are 



(51) 



2v 





1 


2 


3 


4 


5 


6 


7 


T 


3.06362 


2.61931 


2.34104 


2.14351 


1.99308 


1.87315 


1.77429 


1.69088 


2u 


8 


9 


10 


11 


12 


13 


14 


15 


T 


1.61924 


1.55650 


1.50123 


1.45180 


1.40735 


1.36697 


1.33003 


1.2963 


2u 


16 


17 


18 


19 


20 


40 


200 


2000 


T 


1.2650 


1.23616 


1.20927 


1.18411 


1.16058 


0.87348 


0.4229 


0.13888 



To study T,^ for z/ > 10 define 



1 



Pi 



Proposition 9.1. The sequence Ty is strictly decreasing to 0. For z/ > 10 we have 

27r „ 2tt 



(52) 



<Ty < 



- iPty 



Remark 9.2. The zeros (and hence the eigenvalues = jl) are rather well-known, 
namely 



^ 20 



3 2 3/7T,_i/3 nn«i,,-l ^ „■ ^ ,,1/3 , 3 o 



^1 1/3 

^ 20 



for all V G |N with v > 10. Here, ai ~ —2.33811 is the first negative zero of the Airy 
function Ai(x). Therefore, 

(53) u + a v^l^ + h u~^'^ - cu^^ < < u + a u^'^ + h u'^'^ 

with positive constants a ^ 1.8557, b ~ 1.0331, c < For we obtain the estimate 

(54) z/2 + 2a + (26 + a^) v^l^ + 2ah + z/-'/^ _ < A, < 

+ 2a + (26 + a^) v^l^ + 2a6 + 6^ //"^/^ 
where C(z/) = c (2 + 2a + 26 z/"''/^ + c z/~^) is strictly decreasing, and C(9) < 1/5. O 
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We start with proving the estimate fl52l) . which by fl50|l is equivalent to 
(55) < pu < p+. 

Recall that 

(56) 



h[s) = : — = 2u 



Since Ju^i{ju-i) = we have h^jy^i) = 2u. This and h\s) > on (0,j^) show that 

ju-l < Pu< ju- 

In order to improve these bounds on we need to better understand h on the interval 
'■= [ju~i,ju)- The identities P6|) and fl56|) show that 



(57) 

In particular, 
(58) 

Moreover, using ( !57|) 
h"(s) 



h'(s) = s + -(h-2u). 

s 



h'{p. 



2z/ + 1 



h , . h's — h ^ , 
-h' + — (/i-2z/) 



l + h+{h-2i^) 



h^-h 



+ 1 + A (/, _ 2z.) ) . 



It follows that h" > on I,^. Therefore, the straight line of slope h'{j^_i) passing through 
(j,y_i, 2i/) reaches the height 2z/ + 1 on the left of the graph of h, while the straight line of 
slope h'{pv) passing through {ju-i, 2v) reaches the height 2z/ + 1 on the right of the graph 
of /i, cf. the figure below. 

h{s) 



2v 



2v 




Ju-l 



> Pu Pf 



Pu 
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Together with (|58|1 we conclude that 



1 ■ 1 + 

P. + ^ 3.-1 

Using now that < < ju-i + and that > z/ + ^ for all > 10 by (l53i) . we 
find 

2z/+l 2z/ + 2 

Pi/ H < ju-i H ^ < ju-l + 2, 

and hence (!55l) follows. 

It has been shown in [20] that < whence converges to 0. This also follows 
from fl52l) and 



(59) K-iptr = jI-jI-1-2--^ = 2i. + 0{u'/') 

Ju-l 

where for the last identity we used (154|) . Note that (!59|) and Aj, — (p,;)^ = 2z/ + 0(z/^/^) 
imply that 

2n 



We finally show that the sequence T^, is strictly decreasing. In view of the Table fl5T|) we 
can assume that z/ > 10. By fl52|) we need to show that for each such z/, 

{p~Y < A,+ i - {pt+i)\ 

i.e., 

(60) A,+ i - A, > A,_i - A,_i + (^2 - 2-^^) 



The first bracket on the RHS is equal to 

4 m 4 1 

< < 



Ju-l + 2 v + 2 - 3' 
and the second bracket is less than It therefore suffices to show that 

(61) K^^,-K > Vi-A.-i + i + ^. 

The function z^ i— )■ z/° is convex for a = | and a = — |, but concave for a = |. At z^ = 10 
we have 



2b) ((z. + i)^/^ - u'/') > {a' + 2b) ((z. - i)^/^ - (z. - 1)^/^) - 1. 



Furthermore 
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and C{u - 1) < C(9) < | for z/ > 10. Since | + ^ + ^ + 2| < 1, the estimate §M> now 
implies that ( 16T]) holds true. □ 

Remark 9.3. It is known that the function z/ H- Aj^ is strictly convex on (0, cxd), see [6]. 
In particular, 

A,^+i - A,^ > A^_i - A^_i. 
This is not quite enough to prove inequality (!60l) . 
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